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Abstract 

Motivated by the recent measurements of Kondo resistivity in thin films 

and wires, where the Kondo amphtude is suppressed for thinner samples, 

the surface anisotropy for magnetic impurities is studied. That anisotropy is 

developed in those cases where in addition to the exchange interaction with the 

impurity there is strong spin-orbit interaction for conduction electrons around 

the impurity in the ballistic region. The asymmetry in the neighborhood of 

the magnetic impurity exhibits the anisotropy axis n which, in the case of 

a plane surface, is perpendicular to the surface. The anisotropy energy is 

/S.E = if(i(nS)^ for spin S, and the anisotropy constant is inversionally 

proportional to distance d measured from the surface and > 0. Thus at 

low temperature the spin is frozen in a singlet or doublet of lowest energy. 

The influence of that anisotropy on the electrical resistivity is the subject of 

the following paper (part II). 

PACS numbers: 72.15.Q, 73.50.M, 71.70.E 
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I. INTRODUCTION 



In the last few years the experimental study of dilute magnetic alloys with non-magnetic 
host of reduced dimensions has attracted considerable interest The subjects of most 

of the experimental works are the size dependence of the Kondo effect, thus to deter- 
mine whether the Kondo temperature and amplitude depend on the film thickness or the 
diameter of the wire or not. A very recent paper of N. Giordano |^ has indicated that the 
magnetoresistance above the Kondo temperature depends also on the thickness of the film. 

The theoretical motivation of these experiments has been the concept of spin compen- 
sated Kondo state. Considering the Kondo effect the ground state is a singlet where the spin 
of the magnetic impurity is screened by the spin polarization of the conduction electrons, 
which is known as the screening or compensation cloud [^]. The theoretical studies of the 
Kondo effect suggest, that the size of that screening cloud, ^ is of the order of hvp/^kBTx) 
where vp is the Fermi velocity and is the Kondo temperature That length scale is 
especially large for alloys with small Kondo temperature thus for Au(Fe) it can be in range 
of ^ = lO^A with Tk ~ 1 K. In the case of wires or films it is easy to prepare such samples 
where the size at least in one direction is smaller than that Kondo coherence length ^. The 
question has been raised concerning those experiments where the size is smaller, whether the 
coherence length does prevent the formation of the spin-compensated ground state or not 
imj^. Even if that argument looks very challenging, the theoretical base for that argument 
is very weak, as the magnetic impurity experiences the conduction electron density only 
at the site of the impurity. At zero temperature in case S = 1/2 the polarization cloud 
must contain one electron with spin antiparallel to the local spin. The decay of the cloud is 
determined by the correlation function (Scr(r)) where S stands for the impurity spin located 
at r = and cr(r) is the spin polarization of the conduction electron. In dimension d = 3 
the correlation (S<T(r)) decays like r^^, but in reduced dimensions the decay is weaker, for 
= 2 (c? = 1) it is like r^^ (r°), the Kondo coherence length is, however, not affected 0. 
Beyond the distance ^ the decay is exponential like exp{— r/,^}. Thus in ci = 3 the shape 
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of the cloud is sphere, in d = 2 pancake, and in d = 1 cigar hke (see Fig. |l]). Following 
that argument only the level spacing AE can hinder the formation of the compensated 
groundstate if AE > ksTK- That can happen only for grains of size about lOOA, but that 
is out of question for films and wires if the electrons are not localized |0. Thus accepting 
the existence of that size dependence an other explanation is required, but there are also 
experiments where the existence of the size dependence is questioned @. 

The infiuence of non-magnetic impurities on the formation of the Kondo resonance has 
been investigated for thin film using the theory of weak localization 0. In contrast to that, 
the present work deals with the ballistic region. 

Recently it has been suggested by B. L. Gyorffy and the authors of the present paper 
PI that for the magnetic impurity interacting with the conduction electrons by the effective 
exchange interaction a surface magnetic anisotropy can be the result of spin-orbit scattering 
of the conduction electrons on the non-magnetic host. The magnetic anisotropy energy is 
given for a single impurity by the formula 

= MnSy, (1) 

where (nS) is the spin component of the impurity spin in the direction parallel to the normal 
vector n of the surface and the amplitude Kd > is inversely proportional to the distance 
of the impurity d, measured from the surface. 

The magnetic anisotropics caused by the relativistic corrections occuring in the Dirac 
equation, as the dipole-dipole and the spin-orbit interactions can reflect the geometry of 
the sample. For example in magnets they are responsible for the easy axis magnetization 
where the dipole-dipole term is dominating 0. In the case of superimposed magnetic and 
non-magnetic layers a magnetic anisotropy is developed from flrst principles which is formed 
as a result of competition between those two interactions pO|JTl[| . 

As far as it is known by the present authors until recently the possibility has not been 
explored that the spin-orbit interaction between the non-magnetic host atoms and the con- 
duction electrons can produce a magnetic anisotropy for the magnetic impurities by the 
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exchange interaction between the impurity spin and the conduction electrons. Such an 
anisotropy cannot develop for the impurity in the bulk, but that can exist in host limited 
in space. Thus, that anisotropy reflects the geometry of the sample and the position of the 
impurity in that. 

The present paper is devoted to calculate that anisotropy in the second order both in the 
exchange interaction between the impurity and the electrons and in the spin-orbit interaction 
between the electrons and the host atoms. The mean fleld calculation does not lead to such 
terms. 

The present paper is organized as follows. In Sec. |I| the model is described where it is 
assumed that the spin-orbit interaction takes place on the localized e.g. d-levels of the host 
atoms and that localized orbital is hybridized with the conduction electrons following the 



idea of the Anderson model |jT2|. In Sec. |T| the electron Green's function is calculated in 
the flrst order of the spin-orbit interaction. Sec. |I^ is devoted to calculate the impurity spin 
self-energy in which the spin anisotropy given by Eq. (|1]) shows up. The expression for the 
anisotropy constant is developed in Sec. |V| [|T^. The conclusion is presented in Sec. |V|. In 
Appendix ^ the role of the time reversal symmetry is explored. The complicated integrals 
appearing in the flnal expression of the anisotropy constant are developed by analytical 



calculations in Appendix y. The following paper. Part II [|14| deals with the calculation of 



the amplitude of the Kondo resistivity anisotropy as a function of the fllm thickness [13 



The study of the magnetoresistance is left for a further publication ||T5[. The comparison 



with the experiments is contained by Part II [O]. 



II. THE MODEL 

For the sake of simplicity we consider an inflnite half space where the host atoms with 
spin-orbit interaction are homogeneously dispersed (no crystal structure effect), and the 
magnetic impurity is placed in a distance d from the surface (see Fig. For further 
simpliflcation the shape of the sample is taken into account only in the positions of the 
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host atoms representing the spin-orbit interaction, while the free electron like conduction 
electrons move in the unlimited whole space. 

The interaction between the conduction electrons and the magnetic impurity is described 
by the simplest realistic Hamiltonian with orbital quantum numbers. Therefore 5 = 5/2 
is chosen, because in this case the Hamiltonian is diagonal in the orbital quantum number 
m according to the Hund's rule. For other values of S the Hamiltonian consists of several 



complicated terms [[16i[r7| . After these considerations we can write the Hamiltonian as 

Ho = ^^kdl^^akma 
k,m,a- 

+ J X! S {al,^^(T„„:ak'm'a') ^mm' (2) 
k,k' ,m,m' 

where ali^^ (o-kima) creates (annihilates) an electron with momentum k, angular momentum 
/, m and spin a, J is the effective Kondo coupling, cr stands for the Pauli matrices and the 
origin is placed at the impurity site. Keeping only the / = 2 channels the index / is dropped. 

In order to make transparent that the spin-orbit interaction is mainly in the d-channel 
and due to the Coulomb potential of the nuclei we introduce a simple model where the spin- 
orbit interaction takes place on the d- levels of the host, which hybridizes with the conduction 
electrons. These host atom orbitals are labeled by n referring to the position R„ and also 
by the quantum numbers /, m, a (e.g. / = 2 for Cu and Au host and the index / is dropped 
again). The Hamiltonian of these extra orbitals is 

Hi = Bo EbXl + A E HL\m'){a\cT\a')btl^b^i, 



nmm 



+ E (Vkmmi'Rn)b^;:l^akm'a + h.c)j (3) 
nkmm'a 

where b^^l^ (bml) creates (annihilates) the host atom orbital at site n with wave function 
^l=2,m- H;mm'(Rn) is the Audcrsou' hybridization matrix element [^], which depends on 
R„ since spherical wave representation with origin at the magnetic impurity is used, A is 
the strength of the spin-orbit coupling, and L is the orbital momentum at site n. As the 
spin-orbit interaction is weak, therefore, its effect will be considered as a perturbation. 



The Rn dependence in the Anderson' hybridization matrix element can be evaluated as 

nnw(Rn) = {^t2\Heff\^lvn.'), (4) 

where H^ff is the effective hybridization Hamiltonian between the host atom orbital {^^i^) 
and the conduction electron states with origin at the magnetic impurity After 
inserting a complete orthonormal set of free spherical waves with origin at the host atom at 

-Rn 



d^' E I^S"V')(^K'V"I = 1' (5) 
,// // 

and taking into account the usual assumption 

SO that the hybridization matrix element is diagonal in quantum numbers /, m, shows slow 
/c-dependence and it is the same for each host atom we got for the hybridization matrix 
element 

oo 

^fc/zw'(R„,) = Vij rfA;"($K„|$°,„,). (7) 



Thus to get the hybridization matrix element one has to calculate the overlap between 
spherical waves with different origin, 

(*i"Ll*°^'™') = / dv%..,^{v)^u'm'{r + RO (8) 

where ^kim{'^) = \J'^ji{k^)Yi^{~) are free spherical waves, ji{kr) are spherical Bessel func- 
tions, Y{^{r) are spherical harmonics, and r = |r| [|l^]. This can be simplified by using a 
local coordinate system for each host atom where the zi^n) axis is directed parallel to R„. In 
that system m is conserved and 

Vkll'mm'iRn) = ViSmm'Vkll'miRn) i (9) 

where 



VkU'm{Rn) — Sr 



(2/+ -m)! 



'7r\ 2(/ + m)! \ 2(/' + m)! 



(2/' + l)(/'-m)! 



kRr, 







and 



(10) 



Si = j dxxji{x). (11) 



In Eq. (p!0|) the y = r / Rn, y' = |r + Rn|/-R„, in Eq. ( pA]) the x = k"r new integration variables 
have been introduced, are Legendre polynomials, and Rn = |Rn|- The occuring integrals 
could be evaluated analytically. 

For Z = Z' = 2 5*2 = 2, thus by introducing the notation V = V2 



(12) 



where the Vkm{Rn) matrix elements are symmetric for ±m and show different power behav- 
iors in Rn dX k = kp'- 

(D \ _ .^/^ sin(fcf^i?„) 3 cos{kpRn) _ Usin^kpRn) _ 27 cos{kpRn) 27sm{kFRn) \ 
VkM n)- y 2fc^i?„ + {kM' {kpRnf {kpR^Y + {kpRnf )' 



,^ , / cos{kpRn) bsmikpRn) 12cos(kpRn) 12sm(kpRr, 



(kpRn)^ {kpR. 



{kpRr, 



(D \ - . t-/ sm{kpRn) _ 3 cos{kpRn) 3sm{kpRn) \ 
Vk,2[Kn)-i^^ (A;^i?.)^ (kFRnr {kpRnf J' 



(13a) 



(13b) 



(13c) 



III. ELECTRON PROPAGATOR IN FIRST ORDER OF SPIN-ORBIT 

COUPLING 



The electron propagator leaving and arriving at the impurity was calculated in first order 
of spin-orbit coupling according to the diagram shown in Fig. |^. 
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In the local system the conduction electron propagator has the following matrix form in 
first order of spin-orbit coupling 

kma,k'm'u'\1'^n) = : r VV mm' aa' [-n-n ) [i-^) 

lUJn - Sk ^ lUJ - Ek lUJ - Sk' 

where 

Wmm'aa'{Rn) = Vk^mmG d{uj) \{m\L\m') {a\(T\a') G diuj)V^^^,^, . (15) 

The scatterings on several host atoms give higher order corrections in A. 

In Eq. (|15D the ^-dependence was replaced hj kp and Gd{uj) denotes the electron prop- 
agator for the d-levels of the host atom. Gd is given by the spectral function pd as 

Gdiu;) = I P^dJ (16) 
J to — to 



where 

1 A 



(17) 



and A = tiV'^Pq is the width of the d-levels due to the hybridization po is the density 



of states of the conduction electrons for one spin direction. For uj <ti maxje^, A} Gd can be 
replaced by a constant I/eq. 
Thus using Eq. (|l|) 

\T/2 / s 

W^m'aa' (Rn) = B+ a' + a+ + B'a' (18) 

^0 \ / mm' aa' 

where and are 5x5 matrices in the quantum number m, having the form 



= \/(3 + m')(2 - m')vkpmVkFm'8m,m'+i, (19a) 



-Bmm' = \/(3 - m')(2 + m')VkpmVkFm'Sm,m'-l, (I9b) 
^mm' = mVkpmVkpm'Sm^m'- (I9c) 

These matrices could be introduced phenomenologically also. 
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By rotating back the local system to the frame of the sample where the z axis is perpen- 
dicular to the surface the electron propagator can be calculated. These rotations were done 
in the standard way by using the formula 

Wmm'o-o-' {Rn, On, fn) = 

R^Xi^n, 0)4V'^((/^n, ^n, 0)W^^^,^^, (0, -^„, -ipn)R^^C) {Q, -^n, -^n), (20) 

where (i?„, 6'„, ipn) are the polar coordinates of the host atom labeled by n in the system 
of the sample and R^'^\ R^^'^^ are the rotation matrices with angular momentum J = 2 and 
1/2, respectively |T^. In a case when rotation symmetry around the z axis of the system of 
the sample is obeyed, as in the model described in Section |I|, the electron propagator does 
not depend on the azimuthal angle ipn and thus it can be written as 

^mm'(Ta'{Rn,dn) = Sm+a,m' +a' dmm{(^n)d^aI'^\On)Wfnm' aa' (21) 

where the Wigner-formula for rotation matrices [19| was used. 

The time reversal symmetry gives restrictions for the electron propagator (see Ap- 
pendix ^ which provide a check of calculations. In the calculation the angular dependences 
are very important because in the case of s-wave scattering the spin-orbit interaction cannot 



influence the dynamics of the impurity spin . 



IV. SELF-ENERGY CORRECTIONS FOR THE IMPURITY SPIN 



The self-energy was calculated by using Abrikosov's pseudofermion representation [21] for 
the impurity spin and Matsubara's diagram technique applied for the exchange interaction 
with coupling strength J given by Eq. (^. It can be shown that the Hartree-Fock diagram 
gives no contribution. 

The diagrams for the self-energy of the impurity spin which contain the electron prop- 
agator calculated in Section |T| are shown in Fig. |^. The spin factors of these diagrams 
are 
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J SMM''''(T3(Ti'^m3miW^mim2(Ti(T2(-Rn; ^n)SM'A/'''(T2(73^m2r7i3 (22a) 
j'^^MM''^azu-i_^mzmiWm-im2cn(T2i.-^n) ^n) Sm'M''' 0-20-4 '^m2m4 ^m4rra3(T4o-3 (-Rn' ) ^n') (22b) 
•j'^^MM'^cTiai^m4mi^mim2CTicT2{Rnj ^n) W^?Tt2 013 0203 (-^n' ; ^n')^Af'Mf 0-304 '5m3m4 (22c) 

for a), b) and c), respectively. In the spin factor of diagram a) the trace of W disappears. 
This easily can be seen from the form of W in the local coordinate system because 
and are traceless and trace is invariant under rotation. The spin factor of diagram c) is 
proportional to 5*^ thus it does not give contribution to the anisotropy constant. The spin 
factor of the remaining diagram b) is 

S'^Fi{Rn, 9n, Rn', On') + SlF2{Rn, dn, Rn' , Gn'), (23) 



where 



Fl{Rn,On,Rn',(^n') — 2J^/ll22, (24a) 
F2{Rn, On, Rn', On') = 2J^(/llll — /l212 — /1122), (24b) 



with 



/0-10-20304 ^ ^ ^^mm'oi02 (-^nj ^n) ^^m'mo403 (-Rji' ; ) ■ (^^) 
mm' 

After calculation of the remaining part of the diagram Fig. ^ (b), its total contribution 

is 

^PoDf (^^) {S^FiiRn, On, Rn', On') + SlF2{Rn, On, Rn' , On')) (26) 

where po is the density of the states of the conduction electrons for one spin direction and 
D its band width. The function / gives the analytical part of the diagram which is given in 
Table |. 

As this diagram contains two host atoms, averages have to be taken over n and n' which 
will be performed in the next Section. 
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V. THE ANISOTROPY CONSTANT 



As it was shown in Section the leading contribution in spin-orbit couphng to the 
anisotropy constant (see Eq. (P) comes from a second order diagram, namely from Fig. ^ 
(b) and it is 

^ptDf (^) F^iRn, On, Rn', On')- (27) 

As that diagram contains two host atoms with indices n and n', the summation over 
those must be carried out. According to our simple model this gives for the anisotropy 
factor 

K=^J d'Rn J (i'Rn'AplDf F2{Rn, ^n, Rn' , On'), (28) 

where is the size of the volume per host atom. The integrations were calculated by 
considering first the shells with constant i?„ and Rn' (see Fig. |^) and integrating with 
respect to the angles. The integration with respect to (fn and (fn' was trivial according to 
the conservation of the z component of the angular momentum which was used from the 
beginning. 

If e.g. Rn > d then the presence of the surface appears as a limit in the 6n integra- 
tion, more precisely we have to integrate from On,min = arccos((i/-R„) to tt. The integrals 
were calculated in a way in which the integration regime was divided into four parts where 

Rn' < d; Rn > d, Rn' < d; Rn < d, Rn' > d and Rn, Rn' > d, respectively. The integra- 
tions with respect to 9n and 6'„/ was simple and made the contribution of the first part to 
be zero. The others give 

_^ oo d ^ oo oo 

I ^"^""^"/ ^^ri'Rn'Jl{Rn,Rn') + —^ J dRnRn J dRn' R?n' -hiRn, Rn')^ , (29) 

d ro d d 

where ro is a short distance cutoff in range of the atomic radius, and 

TT TT 

Ji(i?„, Rn') = 2 j dOn sin On J dO 
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J2{Rn,Rn')= j d9nsm9n J dOn' Slli 9^' F {R^, Rn' , On, 9^') ■ (30b) 

These remaining integrations with respect to i?„ and were estimated in the leading order 



in l/{kpd) (see Appendix iBj). It turned out (see Eq. ( P20|) and (P22 )) that the dominant 
contribution arises from the integral of Ji (/?„,/?„/) where Rn > d > Rn' or the opposite, 
and the contribution comes from the lower limits of the integral (see Fig. |^), namely Rn = d 
and Rn' = tq. Thus for kpd':^ \ 

,,.,«.„,„,^;(^)_J_f^>0 (3X) 



where P{kprQ) is a numerical factor depending strongly on tq (see Eq. ( |B21D ) and it is 
positive at least for kpTQ > 0.1. 



VI. CONCLUSIONS 

It is shown in the present paper that for a magnetic impurity embodied into an infinite 
electron gas a magnetic anisotropy given by Eq. (0) is developed if the impurity is sur- 
rounded by atoms (e.g. Au, Fe) with large spin-orbit interaction in an asymmetrical way. 
The condition for formation of that anisotropy is that electrons scattered by the magnetic 
impurities are in angular momentum channels different from zero (/ 7^ 0). In the other case 
(/ = 0), the impurity experiences the host atoms in the same distance from the impurity in 
an identical way, thus the shape of the sample does not play any direct role and, therefore, 
no anisotropy axis can be exhibited. 

We considered an infinite half-space with homogeneously dispersed host atoms with spin- 
orbit interaction, in which an impurity is placed in a distance d from the surface of that 
half-space. The shape of the sample was taken into consideration only in the position of the 
host atoms, so the conduction electrons were assumed to move in the whole space. In the 
calculation no randomness was taken into account, therefore, it is valid only in the ballistic 
region. 
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To describe the interaction between the conduction electrons and the magnetic impurity 
we used the simplest realistic Hamiltonian with orbital quantum numbers. 

For the spin-orbit interaction taking place on the d-levels of the host an Anderson like 
model is developed with the spin-orbit interaction on the atomic level of strength A 



and the hybridization matrix element V. In this way for the effective spin-orbit interaction 
between the conduction electron and the host atom an oversimplified model is obtained (see 
Eq. ([TsD). The exchange interaction and the spin-orbit interaction was assumed to be weak, 
thus perturbation theory was applied (see Sec. H). 

First we calculated the electron propagator in first order of spin-orbit coupling (see 
Section pp . The angular dependence was very important because keeping only the s-wave 
scattering the spin-orbit interaction cannot influence the impurity spin dynamics [pO|j . 

Then the self-energy corrections for the impurity spin were calculated by using 
Abrikosov's pseudofermion representation for the impurity spin and final temperature 
Green's function technique for the exchange interaction (see Sec. [TVI) . It turned out that the 
first correction to the anisotropy constant defined in Eq. (|I|) comes from a second order dia- 
gram in spin-orbit interaction (see Fig. ^ (b)), thus an average had to be performed over the 
two host atoms (see Sec. |V| and Fig. ^ for the relevant regions of integrations). The result of 
this averaging was performed in the leading order in l/{kpd) in Appendix 0, the anisotropy 
constant obtained behaves like ~ l/lkpd) in the leading order, and oscillations occur 
only in the next order. 

This behavior turned out to be independent of the actual nonzero values of the angular 
momenta when the calculation was repeated for different angular momenta of the magnetic 
impurity (/') and also of the dominant spin-orbit scattering channel at the host atoms (/). 
Furthermore, in all of the cases K^i > 0. 

To estimate the order of magnitude of the anisotropy factor given by Eq. (|3lD, we consid- 
ered the parameters as Jpo ~ 0.1 (in the case of relevant Kondo temperature Tk ~ 0.1 — 1 
K), A ~ 1 eV, A ~ L) ~ 5 eV, £o ~ 2.5 eV, F{uj/D) ~ 2 (see Table kpa 3 and 



P{kprQ) ~ 10 — 700, kpTQ ~ 0.3 — 1.5 |£2l- Thus, the final estimation for its order of 
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magnitude is 

eY <K < -4t- eV. (32) 
In the present theory for the anisotropy the following approximations are made: 

(i) The electrons form an infinite sea and only the distribution of the spin-orbit scatterers 

reflects the " shape" of the sample. In a real sample conduction electrons are confined 
into the sample and they can scattered by the surface. In a mesoscopic sample the 
surface scattering is rather incoherent because of the absence of smooth surface, there- 
fore, we expect that the qualitative results are not sensitive on the particular model 
considered. 

(ii) It is assumed that the electrons scattered by the magnetic impurities do not change their 

azimuthal quantum number m. That assumption is valid only for perfectly developed 
S = {21 + l)/2 spin {S = 5/2 for / = 2) p|,|l3 (see Sec. |I| and the Hamihonian 



given by Eq. (Q)). If S* 7^ (2/ + l)/2 then the Hund's rule does not ensure that 
conservation and the Hamiltonian must contain several terms. We believe, however, 
that the quantitative result including the relative distances between levels can be 
affected by that generalization, but the concept of surface anisotropy remains valid. 

(iii) The electrons are treated like free electrons, thus their elastic mean free path l^i = 00. 
In the reality the elastic mean free path is finite and the Green's function connecting 
the impurity and the spin-orbit scatterers contains an exponential decay. That decay 
factor ensures that the anisotropy is influenced only by those spin-orbit scatterers 
which are inside of the region of elastic mean free path. 

The geometry treated in the paper is the most simple example. The situation is somewhat 
more complicated e.g. in case of such a thin wire where even the middle of wire is affected 
by the anisotropy. Then nearby the surface the anisotropy axis is parallel to the normal 
direction of the surface, in the middle of the sample the spin direction in the ground state 
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for S = 2 must be, however, perpendicular to all of the surface elements, thus it must lie 
along the axis of the wire. That corresponds to an anisotropy in the direction of wire but 
with a negative coefficient. 

Considering the experimental verification of the surface anisotropy, there are no direct 
evidences. Recently Giordano ^ has performed an experiment which proved difficult to 
explain with previously existing theories and he proposed the presented anisotropy as a 
possible proper theoretical explanation. In that work the magnetoresistance of a thin 
film is studied well above the Kondo temperature as a function of the external magnetic 
field. It was found that the thin film samples needed larger magnetic field to saturate the 
impurity contribution to the resistivity. Considering the surface anisotropy in the presence 
of the field B, the Hamiltonian of the magnetic moment is 

H = K,iSn' + BgBS\ (33) 

where the field is perpendicular to the film. The levels e.g. for 5* = 2 as the function of the 
magnetic field are shown in Fig. |^ without and with surface anisotropy. 

It is clearly shown in Fig. |^ that in the presence of the anisotropy, because of the level 
crossing, a larger field B is required to separate the lowest energy state from the other levels 
in order to saturate the magnetoresistance. The detailed theory will be published elsewhere 

The application of the present theory for resistivity of samples where the anisotropy 
affected regions of impurities are not negligible compared to the bulk, are the subject of 
the following paper [|14| (Part II). In case of S" = 2 at the surface the Kondo effect cannot 
develop as the spins are frozen in the state = 0. 

The theory could be developed further in different directions to include the surface scat- 
tering, to determine in the framework of a realistic atomic calculation, to consider dif- 
ferent geometries, and to take into account the elastic mean free path l^i in an explicit 
form. 

The strongest ambiguity in the calculation is the short range cutoff tq appearing in Eq. 
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(iH). 

Finally it is important to emphasize that the present calculation is beyond the Hartree- 
Fock approximation (see Sec. |IV|), thus that anisotropy should not be obtained by band 
structure calculation in agreement with the present results obtained by L. Szunyogh et al. 
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APPENDIX A: 

In this section we derive the restrictions on the electron propagator calculated in Sec- 
tion ^ due to the time-reversal symmetry. As in we calculate 



^ a 

= I E e"^""" (^a|Q„.(r, t)\^>,) {^M'^'.' {r, 0) I VI/,) 

a,b 

= |Ee-"''"^^^'''""''^^^<^|cw(r,0)|vI/,)(vl/,|cLv'(^,0)|vI/.), (Al) 

a, 6 

where cima{r,t) {clj^^{r,t)) annihilates (creates) a conduction electron of spin a and orbital 
momentum l,m at position r and time t. When the time-reversal symmetry is obeyed 

{ciUr, t)4^'Ar, 0))= {K+CiUr, t)KK+4^,Ar, 0)K) 

= ^ E e-^^"+^(^'"-^^^v[;Jir+Q™.(r, 0)K\^,){^,\K-^cl^,Ar, 0)K\^a) 



z 

a. 
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lf"^'''^'^'^'{^a\cv,-n.'Mr.mb). (A2) 



where K is the time-reversal operator. 

In comparison with Eq. (^l|) we obtain the relation 



(cw(r,t)cL..,(r,0)) = (-l)i+'^+'^'+'+''+-+™'(Q,_,_.,(r,^ (A3) 

Applying the same procedure to (cJ/^/^/(r, 0)Qm.<T('", ^)) the obtained relation is 

(cL'.'(r, 0)Q_(r, t)) = (-l)^+'^+'^'+'+''+™+-'(4_,_(r, 0)c,,_, _..(r, t)). (A4) 

Thus in the case of s-wave scattering (/ = /' = m = m' = 0) the electron propagator is 
diagonal in spin space in agreement with [^. In our case (/ = /' = 2) the restrictions for 



the electron propagator given by the time- reversal symmetry (see Eq. (|A3|) and Eq. (|A^)) 



are 



Wmm'aa'{Rn, ^n, V^n) — ( — 1)"'^+'^+'^ +m+m ]/[/_^, _^(^Jl^^ Q^^ i_p^^_^ 

which served a good check for the calculation. 

APPENDIX B: 

Here we estimate the integrals 

Il = — j dRnRn J dRn'Rn,Jl{Rn, Rn') (Bla) 
d ro 

h = — J dRnRn J dRn'Rn'J2{Rn, Rn') (Bib) 
d d 

appearing in Eq. (|29|) , in leading order in l/lkpd) for kpd ^ 1. This calculation is very long 
for I2, but similar to Ji, thus we present here the estimation only for Ji, but at the end we 
give the form for I2, too. 

After the integration with respect to On and 6n' 
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■ {3VkpoiRn)VkFliRn) - vlpliRn) + 2Vkpl{Rn)VkF2{.Rn) " '^vl^2i.Rn 

■ {QVkMRn'HARn') + vl^liRw) + 4:Vkpl{Rn')VkARn') + 4:vl^2{Rn')) ■ (B2) 

Substituting the Vkpm{Rn) matrix element from Eq. (|13D into the integral, using trigonomet- 
ric identities and introducing the dimensionless integration variables s = kpRni t = kpRn' 
and notations x = kpd, xq = kpVQ, the integral has the form 

16 2,27rV2A 2,900 2 1 



00 



2 2x/225 225 cos 2s 225 sin 2s 18 207 cos 2s 114 sin 2s 

X 

15 327 cos 2s 39 sin 2s 1 3 cos 2s 1 sin 2s 
^ 8s7 8^7 4s6 ^ 4^ ^ 2s5 ^ 8s4 

^ /315 315 cos 2t 315 sin 2t 45 585 cos 2t 165 sin 2t 
\W 2t8 f ^ 2t6 + 2^6 + 

xo 

15 495 cos 2t 63sin2t 1 21cos2t lsin2t\ 
^ 8t4 8t4 4t3 +8t2+ 8t2 + 4t ^ 

The occuring integrals are the type of 

y cis(s^ - x^)—^ (B4a) 



and 



/3 ^ 

ds{s'-x^)^^ (B4b) 



^ 1 

Jds{s'-x')- (B4c) 



^ , sin2t , 
dt^^ (B5a) 



a 



, cos2t , 
dt^^ (B5b) 
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13 

,11111 , , 

dt— = 7-^- B5c 

t" n-1 n-1 ^ ^ 

Let us consider the first two integrals of the second type. After integration by part they are 

m 



and 



where 



^ ,M - (-i)-'p^{|(-i)-(2.):(||^ 
+ g(-i)'(2. + D! Jp,} 



j=0 

22fe-l 



- ^s{2k + 1, 2t) = (-1)^+1 |X:(-iri(2^ + 1)! 



dt 



2j+2 



o2fc 

+ (-1)'7;T7TT« |2«] (B7) 



22fc-l 

+ (-l)^(^^^^[2t], (B8) 



- |;(-1)-(2>)!7|Pt} 



22fc 



oo 

« [t] = -J du^^ (BlOa) 
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Si [t] = du 



smu 



u 



(BlOb) 



are the Cosine and Sine Integral functions [23|. To estimate the integrals in Eq. (P5|) for 
a or (3 ^ 1 the expression of the Cosine and Sine Integral function in terms of auxiliary 
functions was used 103 



Ci [t] = f{t) sint - g{t) cost 



(Blla) 



71 



S^[t] = - ~ fit) cost ~g{t)sm{t), 



(Bllb) 



where 



1 2' 4' fi' 



i=0 



(B12a) 



, , 1 ^ 3! 5! 7! 



1=0 



f2i 



(B12b) 



It can be seen from these asymptotic expansions that the primitive functions in Eq. (B3) 
are e.g. for a 1 



2k-l 



{2k-l)\ 



sin2« - (2i)! 



2a 



i=k 



cos2a ^ ^ (22 + 1)! 



(2a)2 



i=fc-l 



2i 



(B13) 



ifis{2k + l,2a) ~ (-1)* 



)2fe r 



{2ky. 



cos 2a 
2a 



E(-i) 



sin 2a ^ , 



. (20! 
(2a)2^ 



(B14) 



and 



¥.e(2A;,2a)~(-l)^' 



i2fc-l 



(2A;- 1)! 



vr cos 2a 



sin 2a ^ 



2 2a 
(2a) 2^ 



E(-i) 



. (20! 
(2a)^ 



(B15) 
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i2fc 



^,{2k + l,2a) ~ (-1)* 



(2A:)! 



sin 2a 



2a 



E(-i) 



(2a)2* 



cos 2a; ^ 



{2i + l)\ 



(2a)2* 



(B16) 



Thus in our case when a = x ^ 1, [3 = oo oi a = xq, [3 = x ^ 1 the contributions of these 
integrals in leading order in 1/x are 



dt 



sin 2t cos 2x 



2x" 



(B17a) 



(it 



cos 2t sin(2x) 



(B17b) 



dt 



sin 2t 



-iPs[n,XQ) 



(BlTc) 



(it 



cos 2t 



-V2c(n,xo), 



(B17d) 



where 99<j(n, xq) and <y9c(^, a;o) denote the primitive functions of the integrals given in Eq. 

(0), (0), o, (0). 

Turning to the integrals of the first type in Eq. they can be transformed by inte- 
gration by part into 



/3 



dsi^s^ — x^ 



sin2s Is^ — x"^, cos2sJ^ 



P 13 

, , cos2s f , n , n 9, cos 2s 
+ / ds / ds-(s^ - x^) — 

a 



, (B18a) 



ds{s^ — x^ 



cos2s fs — X sin2s1'^ f . sin 2s r , n 



s" 2 



Sill 25 



X 



(B18c) 



.(3-n)s"-3 (l_^)sn-l 

Using the leading order formulas for the integrals of the second type in Eq. ( |B17| ) and 
considering our case {a = x ^ 1, f3 = oo) the integrals above in the leading order in 1/x 



are 



21 



sin 2s sin 2x 



2x' 



i-i ' 



(B19a) 



cos 2s cos 2x 



2x 



n—l ' 



(B19b) 



ds{s^ — x^) 



1 



1 



1 



1 



(n-3) {n-l)Jx'^-^' 
Thus the final estimation for the Ii integral in the leading order in \/x {x = kpd) is 



(B19c) 



/i = 4J^ 



^0 



2\\ 2 



P{xo) P{xo)sm{2x) 6750(1 + cos(2a;)) 



X 



x^ 



x^ 



(B20) 



where 



315 

P(xo) = 54000 (—vPc(8,a;o) 



585 



495 



21 



(Pc{Q,Xo) + —ipc{4:,xo) - —(pc{2,xo) 



63 1 

+ 315(/7s(7, Xq) - 165(/7s(5, Xq) + — <^s(3, Xq) - -<^s(l, Xq) 



45 

2Xn 



54000 



9 5 1 

2^0 8,To 8,To , 
45 9 5 



1 



\ 2 Xq 2 8 Xq 8 Xq 



45 C0S(2 Xq) 



2x1 



81 cos(2a;o) 53 cos(2a;o) cos(2a;o) 45 sin(2a;o) 



\xl 



+ 



2xe 

21sin(2a;o) 5 sin(2a;o)\ 

4 4^i~/ 



8x0 



^6 



The I2 integral in leading order in 1/x {x = kpd) is 

'nV^W" 1 



el 



{kpaf 



6750 cos (2x) ^ 347625 
x^ 64x^ 



(B21) 



(B22) 
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d=3 



FIGURES 




d=2 




d=l 



FIG. 1. The Kondo sreening cloud in different dimensions. 




FIG. 2. The infinite half space of homogeneously dispersed host atoms with spin-orbit interac- 
tion and the magnetic impurity in a distance d measured from the surface. 



1 m o 1 m o 1 m o 1 m o 
• • o • • = • X • 

J V V J 

FIG. 3. The electron propagator leaving and arriving at the impurity. The heavy lines represent 
the localized d-electron propagators, and V and A indicate the hybridization with the localized 
orbital and the spin-orbit interaction, respectively. The indices are according to the local system. 
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n n 




c) 

FIG. 4. The self-energy diagrams for the impurity spin. The double line represents the spin, 
the single one the conduction electrons. The solid circles stand for the exchange interaction and 
the X labelled by n for the effective spin-orbit interaction on the orbital of the host atom at R„. 




FIG. 5. Carrying out the average over the homogeneously dispersed host atoms, the different 
shells with constant radiusis Rn and Rn' are shown which can be restricted by the surface. 




FIG. 6. The two dominating regions contributing to the double integral. Region 2 is formed 
by those shells of smallest radia which are not complete due to the presence of the surface. Region 
1 is around the impurity of smallest radia with short range rutoff tq. 
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a) b) 
FIG. 7. The levels for 5 = 2 as the function of the external magnetic field (a) without and (b) 
with surface anisotropy. 
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TABLES 



X 


ffx) 


n nn 
u.uu 


1 /I Q 1 Q 


U.UD 


i.4yOD 


n 1 n 


i.OD < O 


U. 10 


1 fi47Q 
i .04 < y 


n on 
u.zu 


i. / O ( o 


n OK 
u.zo 


i.oo ( Z 


n Qn 


1 Q/IC/I 


U.oO 


^.UDOo 


n /in 


1 0*^7 
z.iyo / 


U.40 


Z.ooio 


n f^n 

U.OU 


/1 777 


n KK 

U.OO 




0.60 


2.7939 


0.65 


2.9615 


0.70 


3.1334 


0.75 


3.3075 


0.80 


3.4807 


0.85 


3.6486 


0.90 


3.8047 


0.95 


3.9376 



TABLE L The analytical part of diagram Fig. || b) in function of x = 
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